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The energies and widths of DNN quasi-bound states with isospin I = 1/2 are evaluated in 
two methods, the fixed center approximation to the Faddeev equation and the variational method 
approach to the effective one-channel Hamiltonian. The DN interactions are constructed so that 
they dynamically generate the A c (2595) (7 = 0, J^ — l/2~) resonance state. We find that the 
system is bound by about 250 MeV from the DNN threshold, -*/s ~ 3500 MeV. Its width including 
both the mesonic decay and the D absorption, is estimated to be about 20-40 MeV. The 1 = DN 
pair in the DNN system is found to form a cluster that is similar to the A c (2595). 
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I. INTRODUCTION 

The interaction of mesons with nuclei and the prop- 
erty of mesonic bound states are one of the most impor- 
tant topics in the nuclear-hadron physics [iJ-Q- Bound 
states of pions and K~ have been investigated for long 
and have revealed the roles of strong interactions in the 
hadron-nucleus bound states. A step forward in the ex- 
perimental observation of the most deeply bound pio- 
nic states was given using the (c?, 3 He) reaction 0], and 
also, although less clearly, using the coherent radiative 
7r~ capture 0] in Ref. Q. The deeply bound kaon atoms 
had been studied theoretically using the optical poten- 
tials [nJ-Eij]. Because of the large imaginary part, the 
width of the bound states is larger than the energy sepa- 
ration between the levels [EtHI], so that the experimental 
observation is not feasible (see also Refs. (MG3)- 

The simplest of the many-body kaonic nuclear system 
is the KNN, which has had much attention theoretically. 
Because the A(1405) resonance is interpreted as a quasi- 
bound state of the KN system in the nT, continuum [iH - 
Hfl, one expects a quasi-bound KNN system driven by 
the attractive KN interaction in the isospin 1 = chan- 
nel. Various approaches have resulted in a rather general 
consensus that the quasi-bound state is obtained above 
the 7rSJV threshold and the width is larger than the bind- 
ing [26-37]. Thus, the experimental identification of this 
system would be difficult. 

What we report here is the analogous state of the 
KNN, substituting the K by a D meson. The DN inter- 
action in / = is predicted to be attractive in the vector 
meson exchange picture, and thus to dynamically gener- 
ate the J p = 1/2- excited state, A c (2595) ^£3. The 
A c (2595) resonance is rather narrow (T < 1.9 MeV), in 
contrast to the analogous A(1405) with apparent widths 
of the order of 30-60 MeV [13^3 • While the large 
width of the A(1405) is responsible for the large width of 
the KNN state, the analogous state DNN, where the 



A c (2595) plays the role of the A(1405) in the KNN state, 
has much better chances to survive as a long lived and 
observable state. 

The interaction of the D mesons with nuclei has been 
addressed in Refs. [39j, [4l|, |42| and the possibility of mak- 
ing bound atomic states of D mesons in nuclei has been 
considered in Ref. [43| • The interaction of DN is attrac- 
tive both in isospin 1 = and 1 = 1, but much weaker 
for / = 1. Although this leads to a weakly attractive 
D + p interaction, the Coulomb repulsion becomes impor- 
tant for heavier nuclei. As a consequence, the D will 
be only weakly bound in heavy nuclei and the probabil- 
ity to see these bound states is not excessively promising 
43]. However, few-body systems like DNN are less af- 
fected by the Coulomb repulsion particularly for the total 
isospin /tot = 1/2. 

With this in mind we tackle the DNN system from 
two different approaches. The first one is using the fixed 
center approximation (FCA) to the Faddeev equations, 
as done in Refs. [H-[36| for the KNN system. The 
second one is using the variational method as done in 
Refs. [13, H3]. In order to gain confidence that the state 
found is narrow, we have also evaluated the width of the 
state coming from the absorption of the D by a pair of 
nucleons going to the A C N system in the FCA, analogous 
to the absorption of if by a pair of nucleons as consid- 
ered in Refs. [O, H3, H] ■ In the variational approach, 
we extract typical size of the quasi-bound state from the 
obtained wave function. 

The paper is organized as follows. In Sec.|TTl we briefly 
introduce the coupled-channel approach for the DN scat- 
tering and derive the corresponding DN potential. These 
provide the basis of the three-body calculations in later 
sections. The FCA to the Faddeev equations is formu- 
lated in Sec. IIII1 together with the evaluation of the 
two-nucleon absorption. The variational approach to the 
same DNN system is discussed in Sec. IIVI The numer- 
ical results of the three-body calculations are shown in 
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Sec. EJ The discussion for the obtained results are given 
in Sec. I VII The conclusions of this study are drawn in 
the last section. 



II. DN SCATTERING AND INTERACTION 

We consider the two-body DN scattering based on the 
model in Ref. [3i| . This is a coupled-channel approach to 
the s-wave meson-baryon scattering in the vector-meson 
exchange picture. The negative parity A c (2595) reso- 
nance is dynamically generated as a quasi-bound state of 
the DN system in the 1 = channel, just like the A(1405) 
resonance in the strangeness sector [liM35l ]. In Sec. Ill A| 
we derive the DN two-body scattering amplitude which 
will be used in the FCA calculation. An effective single- 
channel potential is constructed so as to reproduce the 
equivalent scattering amplitude in section III Bl This will 
be the basic input in the variational calculation. We work 
in the isospin symmetric limit, which is sufficient for the 
required precision of the present study. 



A. Coupled-channel model for the DN scattering 

We consider seven (eight) coupled channels in the 
isospin I — (I = 1) sector, DN, 7rE c , r/A c , KE C , K"E! C , 
D S A, and rf A c (DN, irA c , 7r£ c , r]E c , KE C , KS' C , D S Y, 
and rj'T, c ). In Ref. [39(, the coupled-channel interaction 
is given by the Weinberg- Tomozawa term 



(W) = - 



kC 



(J) 



4/2 



-(2W-Mi-Mj) 



Mi + Ei Mi 



2M, 



2Mi 



where W is the total energy, / is the meson decay con- 
stant, Mi and Ei are the mass and energy of the baryon in 
channel i, respectively, and Gy is the group theoretical 
coupling strength for isospin /. The reduction factor k is 
introduced to take into account the mass difference of the 
exchanged meson, which we set k = 1 (k = k c = 1/4) for 
the uds (charm) flavor exchange process [39(. The scat- 
tering amplitude tij is obtained from the matrix equation 



(1) 



where the diagonal loop function is given in dimensional 
regularization as 
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FIG. 1. (Color online) S-wave DN scattering amplitude in 
the coupled-channel model |T]) (Left: I — channel, Right: 
1 = 1 channel). Vertical dotted lines represent the threshold 
energies of 7rE c and DN channels. 



where q~i is the magnitude of the three-momentum in the 
center-of-mass frame. We choose the subtraction con- 
stants at fi = 1 GeV as 



a-DN = — 2.056, 



-2.06 (i ^ DN), (2) 



for both the 1 = and 1=1 states, so that the A c (2595) 
resonance is dynamically generated at the observed en- 
ergy. By choosing the isospin symmetric subtraction con- 
stants ©, a resonance state is also generated in I = 1 
at ~ 2760 MeV. The diagonal components of the s-wave 
scattering amplitudes in the DN channel, which are com- 
plex above the 7rY" c (Y c = A c , S c ) threshold, are shown 
in Fig. [1] The resonant nature of the amplitudes can be 
seen in both channels. 

It is worth comparing the 1 = DN-irY> c system with 
the corresponding KN-ttY, system. Both the systems 
have a quasi-bound state. Neglecting the small effect 
of the normalization factor, we can write the coupling 
strength for the DN case as 



2W - M, - Mi 



4/2 



where the channels are assigned as DN (i = 1) and 7rE c 
(i = 2). This is the same form with the KN-ttY, case, 
except for the factor k c = 1/4 in the off-diagonal chan- 
nel. The diagonal interaction is proportional to the me- 
son energy W — Mi, which is reduced to the meson mass 
at threshold. Thus, there are three differences from the 
strangeness sector: 1) heavy mass of D meson, which 
enhances the strength of the DN interaction by the en- 
ergy factor W — Mi, 2) large reduced mass of the sys- 
tem, which suppresses the kinetic energy in the charm 
sector; 3) weak transition coupling DN -4 7tS c , which 
suppresses the decay of the quasi-bound state into the 
7rE c state. These facts explain the reason why the DN 
quasi-bound state is generated with larger binding energy 
and narrower width than those of the KN quasi-bound 
state. In addition, 1) and 2) also enhance the attractive 
interaction in the 1=1 channel. As a consequence, we 
obtain a resonance state also in I = 1 at ~ 2760 MeV, 
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as far as choosing the isospin symmetric subtraction con- 
stants ©. 



B. Effective single-channel DN potential 

Now we construct an effective single-channel potential, 
which will be used in the variational calculation of the 
DNN system. We utilize the method in Ref. (45|, first 
constructing a single-channel framework which is equiva- 
lent to Eq. ([1]) and then translating the result into a local 
and energy-dependent potential in coordinate space. 

The effective interaction v cS is constructed to repro- 
duce the original amplitude in, given by the DN single- 
channel scattering equation (we suppress the isospin in- 
dex in this section) 



(3) 



ft is shown that the v oS is given by the sum of the bare 
interaction in channel 1 (t>ii) and the term with coupled- 
channel effects as [45j 



N N 
cS + ^ v lm g m v ml + ^ Vi m 9mtmi .'// ''/:.• 

m,l=2 



m=2 



(4) 



where = [(v^^)]- 1 - g( N ^)]^ is the (N - 

1) x (jV — 1) matrix of the coupled-channel amplitude 
without the DN channel. In this way, Eq. gives the 
equivalent amplitude with the II component in Eq. (fTj). 
v cS is complex above the irY c threshold, because of the 
imaginary part of the loop function of the ttY c channel in 
Eq. Q. 

We then translate v into the local potential in co- 
ordinate space. Adopting a single gaussian form for the 
spatial distribution, the two-body potential can be writ- 
ten as 



i «n Mn 
v D N{r;W) = , 



uj(W) 

x [v cB (W) + Av(W)] exp[-(r/a 



a) J) 



(5) 



where a s = 0.4 fm is the range parameter of the po- 
tential and uj{W) is the reduced energy of the DN sys- 
tem. The energy-dependent correction term Av(W) is 
introduced to compensate the deviation from the lo- 
cal potential approximation. This complex and energy- 
dependent potential reproduces the scattering amplitude 
tn when the Schrodinger equation with this potential 
is self-consistently solved. The strength of the potential 
vdn{t; W) at r = is shown in Fig. [5] One finds that the 
real part (imaginary part) is larger (smaller) than that 
of the KN potential [45[ , which demonstrate the differ- 
ences of the interaction kernel discussed in the previous 
section. 




2600 2700 2800 
W [MeV] 





-500 
-1000 
-1500 
-2000 
-2500 



• Re v(0,W) 
Im v(0,W) 



2600 2700 2800 
W [MeV] 



FIG. 2. (Color online) Strength of the effective potential 
vdn{v, W) at r = by Eq. (JSJ) (Left: 1 = channel, Right: 
1 = 1 channel). The range parameter is chosen to be a s = 0.4 
fm. Vertical dotted lines represent the threshold energies of 
7rE c and DN channels. 



III. THE FIXED CENTER APPROXIMATION 
FOR THE DNN SYSTEM 

The fixed center approximation (FCA) to the Faddeev 
equations has been used with success in several problems. 
It is advantageous that the two-body absorption process 
of the three-body system can be calculated as discussed 
in Sec. lIII 51 One assumes that a pair of particles remains 
relatively unaffected by the interaction of the third par- 
ticle with this pair. This usually happens when the third 
particle is lighter than the constituents of the pair, and 
also if the cluster is tightly bound. The method has been 
used with success in the study of K~ scattering with 
the deuteron in Refs. (4614491 (see a review in Ref. [5(| 
for comparison with full Faddeev calculations) . More re- 
cently it has been applied to systems of two mesons and 
a baryon in Ref. [51| , where the NKK system is investi- 
gated. The results obtained are in good agreement with 
more accurate results obtained with variational calcula- 
tions in Ref. 15211 , or the Faddeev equations in coupled 
channels [53l. |54| . The puzzle of the A 5 / 2 +(2000) is also 
addressed with this technique, assuming this resonance to 
be mostly built up from irpA in Ref. [551 ] . Closer dynami- 
cally to the problem under consideration is the work [56| , 
where the NDK, KDN and NDD systems are studied 
with this method. 

In the present case, where we want to study the 
DNN system, we have also the precedent of the work of 
Refs. [34l[36| , where the KNN system was studied within 
this approximation and found to provide results in qual- 
itative a gree ment with those of the variational calcula- 
tions [29l |30|| . The condition that the interacting particle 
(D meson) is lighter than those of the two-body cluster 
(nucleon) is not fulfilled in this case. This certainly in- 
troduces larger uncertainties than in other cases studied 
but we still expect that one can get good results at a 
qualitative level. Actually, the real difficulty of the FCA 
occurs when one applies it to studying possible resonant 
three body systems above the threshold of the three par- 
ticles [H3 • In the present case, we look for deeply bound 
states of the DNN system and we are safer. However, in 
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FIG. 3. Diagrammatic illustration of the three-body equa- 
tion ©. 



order to be more certain about the results, we have also 
performed calculations using a variational method. The 
differences found in the two approaches can give us an 
idea of the uncertainties, and the features shared by the 
two approaches can be considered more reliable. 



A. The formalism for the FCA in the DNN system 

In the FCA to the Faddeev equations for the DNN 
three body system, one takes the NN as a cluster and 
D scatters from that cluster. We consider the DNN 
system with total isospin I to t = 1/2 and with the total 
spin- parity J p = CP and J p = IP In this approach, all 
the two-body pairs are in s wave. 

First we make the evaluation for the case of J p = CP, 
which corresponds to the spin (isospin) of the NN pair 
as Snn = (Inn = !■)• To have total isospin I t ot = 1/2, 
the dominant component of the DN system is 7 = 0, 
where the A(2595) resonance appears. 

The T matrix for the three-body DNN scattering is 
labeled by the DN isospins in the entrance channel I 
and the exit channel I', Tjjl We denote the two-body 
(s-wave) DN scattering amplitudes by for 7 = and 
fi 1 ) f or 7 = 1. Then the T matrix satisfies 



Tit, = tW6iT,+tWGiT„G Ti,, r P ex , 



(6) 



which is diagrammatically represented in Fig. B| In 
Eq. ([5]) , Go is the meson exchange propagator [34j, |58| 



Go — 



d 3 q 



q z — m D + te 



(7) 



where FNN(q) is the form factor, representing momen- 
tum distribution of the NN system. P ex is the isospin 
exchange factor, which depends on the total isospin of the 
nucleon, 7jvjv> m the final state, P ex — (— l) /jv «+ 1 = l 
for J = 0, and = —1 for J = 1. 

Here we concentrate on the isospin factors in the DNN 
scattering amplitudes. We define the isospin doublets, 
N = (p,n), D = (D+,-D°) and consider the DNN 
states with the total isospin 7 t ot = 1/2. There are two 
independent states with the total spin J = and J = 1, 
which can be decomposed into the DN isospin eigen- 
states, as 

|^(A r iA r 2)/ NJV =i)j=o = ^|(£Wi) iV 2 ) + ^\(DN-,)xN 2 ), 

1 /3 

|^(^iA r 2)7 NJV =o).7=i = --KDNjoNi) + -y\(DNi)xN 2 ) 



The D exchange matrix is given in terms of the isospin 
recombination factors. 

l(DJVi)oiVa) = i|(DJVa)oJVi) + ^l\(DN 2 )iNi), 

|(UJVi)iJVa) = ^-|(£>JV 2 ) iVi) - -|(DJV a )iJVi). 
Thus the transition matrix G is given by 
/ 1 

G = 




The three-body amplitude T/ j/ is obtained by solving 



Eq. 

T = 



1 - |(t (0) -t^)G Q P ex -t^Gl 

+ It«G t(°)Pe, ^t(°H^G P ex " 
M ^ t (o) t H) GoPex t (D_i t (o) Got (i)p e 



In calculating the T matrix for the scatterings in the J = 
and J = 1 channels, we take the linear combinations, 
with a factor 2 for the choice of the first nucleon, as 



T(J = 0) = 2(f i) 



Too T 01 \ (& 
Tio Tn 

\ I Too T 



r (J = i) = 2(-^) (tZZ (4 



Substituting the T matrix and replacing P ex by +1 for 



J — 0,Inn 
we obtain 



1 scattering and — 1 for J = 1, 7jvjv = 0, 



T(J = 0) = QtW + ItW + 2t^t^G ^j 



T(J = 1) = 



l-^(< (0) -< (1) )G -< (0) i (1) G2 



i t (o) + ^(i )+2t (o ¥ i) Go 



l + i(tW-t( 1 ))G -< (0) i (1) G2 



(8) 



(9) 



These results coincide with those derived in the charge 
basis [33, [36| (see Appendix [B| . 

We can see that Eq. contains the folding of the D 
intermediate propagator with the form factor of the NN 
system. The variable q° in Eq. (J7|) is the energy carried 
by the D, which is given by 



s + m 2 D 



{2M N f 



2V~s 

with y/s for the rest energy of the DNN system. Eq. 
([7]) requires the NN form factor. For 7jvjv = one could 



take the deuteron form factor, but the attraction of the D 
on the nuclcons will make the NN system more compact, 
like in the case of the KNN system. Yet, there are limits 
on how much one can contract this system because of 
the strong NN repulsion at short distances. In order to 
estimate the NN size one can rely upon the results of 
Ref. [3(| in the study of the KNN system, where the 
NN repulsion at short distance was explicitly taken into 
account. In practical terms we use the same expression 
for the form factor as for the deuteron I59l 



F(q) = d 3 P J2 
Jo „•_ -, 



C, 



11 

,2 XI 



c, 



r, (10) 



p 2 + m 2 . 4^ (f— q) 2 + m 
j— i i — i 

but with the parameters wn rescaled such as to give an 
average separation of the nucleons of Rnn — 2 fm [30| . 
They are shown in Fig. [U The validity of this NN form 
factor will be examined by the result of the variational 
calculation, where the average distance of the NN pair 
in the DNN system will be optimized in the three-body 
dynamics. 

We need the argument s\ of the DN amplitude, 
t(y/si). To evaluate it we adopt a common procedure of 
dividing the binding energy into the three particles pro- 
portionally to their masses. The energy of the nucleon 
and the D meson are given by 



E 



N 



M 



N 



2Mm + mo 



E D = m D 



2Mm + m_o ' 



so the total energy of the two-body system can be calcu- 
lated as 



si = (pd +Pn 1 ) 2 = s ( 



M N + m D \ 2 _ ^ 2 
2M N + m D ) Pn °-' 



(11) 



The approximate value of p\ 2 can be obtained by assum- 
ing 



Pn 2 



2M 



D 



N 2 



B N2 = M N - M N 



N 



2Mn + nip 



(12) 



which provides a rough estimate for bound systems with 
the strong interaction. 



B. Evaluation of the D(NN) Absorption 

As we shall see in Sec. |V] we obtain a DNN bound 
system with a very small width. This is related to the 
small width of the A c (2598) state which is generated in 
DN interaction in I = 0. Yet, this calculation only takes 
into account the decay channel DN —> 7rE c for which 
there is little phase space and DN ttA c channel which 
comes from the subdominant DN 1 = 1 component in 
the DNN system. Now we allow the D to be absorbed by 
two nucleons, in analogy to the KNN — > A7V considered 
in Refs. [TEH]. Here the channel will be DNN -4 iVA c 
whose absorption process is shown diagrammatically in 



F(q) 

^^reduced 



1000 
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FIG. 4. Form factor of the deuteron, and the one correspond- 
ing to an iV./V system with a reduced radius from Ref. [30] . 
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FIG. 5. Diagrammatic representation of the D(NN) absorp- 
tion. 



Fig. (other mechanisms and decay channels will be dis- 
cussed in the end of this section). We calculate only the 
first diagram in Fig. [5l The second one gives an identical 
contribution and they sum incoherently: there is no in- 
terference since the A^A C and A C A^ are orthogonal states. 
Hence, the total width will be twice the one obtained 
from just one diagram. 

The S-matrix for the diagram is given by 



d 4 x / d 4 y(-i)t DN ^ DN 



\2uid 1 



-iq(y-x) _ 



q" — m 





xV y aqe iE ^v e ~ iE ^ ^ c (£Wy), 

where V y is the Yukawa vertex. We take the same cou- 
pling as K~p — > A since in the D and A c the c qu ark p lays 
the role of the s quark in the K and A. In Ref. [6fJ, the 
V y is given as 

1 3F + D 
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with D = 0.795, F = 0.465 [H|. We perform the x°, y° 
integrations, and make a change of the spatial variables 



^ t ^ r 

x = R--, y = R+z- 



Then we can write 



ip Nl (x)tpN 2 (y) = 



AP.R,„i-f 



V 



where ip(r) is the wave function of the NN system. N[ 
and A c will be outgoing plane waves. Let us also assume 
that the D is a plane wave with a certain momentum. 
The final formula that we shall use is independent of this 
momentum, as we shall see. Thus, 



1 



v V 



and then, using these new functions, the S- matrix is writ- 
ten as follows 



Q 



1 



V 2 J (2tt)3V2^ 



tDN^DN- 



1 



m n + it 



xV y aq^(q-p Ac + -){2^f5\p % - p f ) (13) 



i T^=^(27T) 4 S 4 (p l -p f ), 



where (p{q) is the Fourier transform of the wave function 
tp{r) normalized to 1, and pi and pf are the initial and 
final momentum, respectively. The NN wave function in 
momentum space is defined as 



and has a maximum value for q — 0. If we take the NN 
system at rest, P = 0, the wave function ip(q) in Eq. (fT5]> 
will peak at q — p Ac — 0. This allows us to approximate 
the D propagator in Eq. (fT3")l as 



(14) 



q 2 — m 2 



D 



(q°) 2 -pI- 



where q° = Ea c — E^ 2 and p^ a w 
X 1 / 2 (Mj lND ,Mj ! ,Ml c )/2M NND . We do not need 
to specify the i_Djv->nA' amplitude since it will be 
accounted for at the end of the formalism. 

Defining of q — p^ c = q ', the square of the total matrix 
element is obtained as follow: 



ITI 



1 



(q ) 2 -p 2 K -m 2 D J 
^2 J q' 2 dq'<p(q')tDN,DN{Vs') 



, (15) 
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FIG. 6. D(NN) absorption. 



With this T matrix we evaluate the cross section for the 
process of Fig. [5] (left) and we obtain 



Cabs 



1 MnnMaJ^pj^ 2 



2tt 



M 2 

1V1 NND 



PD 



It is interesting to relate this cross section to the imagi- 
nary part of the forward D(NN) — > D(NN) amplitude 
from the diagram of Fig. [S] using the optical theorem. 
We find 



Im T D ( NN ) 



Mnn 



Cabs 



1 M A M N 2 

-PK\ T \ ■ 



2n Mnnd 



The next step is to convert the absorption diagram 
of the Fig. [S] (a) into a "many body" diagram of Fig. 
[5] (b) where the nucleon where the D is absorbed, the 
only occupied state of the "many body" system, is con- 
verted into a hole state in the many body terminology 
[62| . Once this is done one observes that if we remove 
the amplitude tr>N hi the expression of T, the expression 
that we obtain for Im Td( NN \ corresponds to the evalua- 
tion of the imaginary part of the two-body loop function 
g of a nucleon and a D meson [Fig. [5] (c)] but with a D 
selfenergy insertion accounting for the (A c Nh) excitation 
of the D meson. We call this Sg. The Feynman rules to 
evaluate Im 8g and Im Td(jvjv) are identical, except that 
toN,DN is removed in the evaluation of Im 5g. Hence we 
obtain 

.1 M Ac M N - 2 
iha.6g = -i— — PaJT\ . 



2tt M 



NND 



with \T\ 2 is given by Eq. (fT5|) removing tr>N,DN- This 
simplifies the expression since 



^ h n dq'(p{q') 



Thus | T | 2 is given by 



m 2 = v v 2 pi 



lim 

r->0 



d 3 q' 
(2tt) 3 

o). 



W') (16) 



[(Q°) 2 -pl-ml} 2 



\m\ 2 



Finally p A 2 accompanying V 2 in the former expression 
requires a small correction. The factor comes from the 




FIG. 7. The meson- baryon loop function goN in the DN 
channel (solid line) and with the effect of the two-body ab- 
sorption ilm Sg added (dashed line). 
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FIG. 8. Modulus of the two-body amplitude DN -> DN 
(solid line) and with the effect of the two-body absorption 
ilm Sg added (dashed line). 



non relativistic aq form of the DNA C vertex. If we take 
instead the relativistic Yukawa vertex of the type 7 At 7 5 , 
then we find the easy prescription to account for the rel- 
ativistic correction, 

KpI -+ v v^r( M » + m aJVa 2 c . 

The next step is to reevaluate the toN.DN amplitude 
used as input in the fixed center formulas. As we men- 
tioned, they were obtained using the method of Ref. [3§| 
with several coupled channels and the formula (JXJ) . We 
redo the evaluation by replacing the loop function in the 
DN channel as 

9dn -> 9dn +Hm8g (17) 

to take into account the D absorption by two nucleons or, 
analogously, the A c Nh excitation of the D meson. When 
doing this, the DN amplitude becomes complex below 
the DN threshold and the narrow A c (2598) resonance 
acquires now a moderate width due to the D absorption 
with a second nucleon. The second process of Fig. [5] 
(right) is accounted for when we consider the three-body 
amplitude T in the FCA formula with the first D scat- 
tering with the second nucleon. 

For the estimation of the width we take the wave func- 
tion 

{\a 2 — q 2 ) 2 + q 2 a 2 ' 

with a ~ 1.7fm , which corresponds to an NN object 
of relative distance 2 fm. 

Let us numerically investigate the effect of the absorp- 
tion using the model described in Sec. Ill Al In Fig. 
we show the meson-baryon loop function gjj /y in the DN 



channel together with the two-body absorption contribu- 
tion to the imaginary part, ilm Sg. We can see that the 
imaginary part of the total g function is no longer zero 
below the DN threshold due to D absorption. In Fig. [5J 
we show the modulus of the two-body amplitude \t\ for 
the DN channel for 1 = using goN and gDN + *I m &9 
of Eq. (fTTj) . As we can see, the inclusion of the absorp- 
tion mechanism induces an increase in the width of the 
peak of A c (2595) in |t| which will have repercussion in 
the width of the DNN system. 

For a narrow resonance, we can approximate the am- 
plitude around the resonance energy by a Breit-Wigner 
form 

t(\/si) ~ p. 

y^sl - Mr + i T 

This leads to the expression of the coupling of the reso- 
nance to the DN scattering state as 

g 2 = l -T\t{M R )\. 

Inspection of Fig. [8j together with the values of T(no 
absorption)= 3 MeV and r (absorption) = 15 MeV, show 
that the value of the coupling g 2 barely changes from 
the introduction of ilm 8g, but of course the reso- 
nance has become wider. Indeed g 2 (no absorption)/ 
5 2 (absorption)~ 6/5. 

The absorption diagram that we have considered is 
not the only one, but it is the most relevant. On the 
same footing we should consider the diagrams where the 
Yukawa coupling produces S c or even E*(2520). The 
analogy with the kaons made before, and the values of 
these couplings that can be seen in Ref. [61] • together 
with the dynamical factor p 3 of the cross section, make 
the contribution of these terms of the order of 5% of the 
A c production and we neglect them. Analogously we can 
also have DN — > 7r(7y)A c (E c ) in the first hadron line of 
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the absorption diagram and exchange a pion or an eta. 
These diagrams are further suppressed because they re- 
quire the exchange of a heavy vector in the DN — > ir(ri)Y c 
amplitude in the extension of the hidden gauge approach 
that we use. They are penalized by the factor k 2 c , which, 
even considering that the pion (eta) propagators have 
bigger strength than the D one, renders these diagrams 
at the level of 10%. 



define the following matrix elements 



171 


— \ V \ 1 ~ J c.m.| * ), 


V{NN) 


=<* |tW|*), 


V(DN) 


= (*|Re VdjvI*), 


T 

-l nuc 


=(*|fjv-T c . m . iiV |* 


Enn 


=T nuc + V(NN), 



IV. VARIATIONAL CALCULATION OF THE 
DNN SYSTEM 



Here we calculate the energy of the DNN system with 
a variational approach formulated for KNN system in 
Refs. [1^, [3(3]. As in the case of the FCA, we consider 
the DNN system with total isospin I = 1/2 and the 
total spin-parity either J p = 0~ or J p = 1~ . The trial 
wave function for the J p = 0~ state is prepared with two 
components: 



7=0 



where J\f° is a normalization constant and C° is a mixing 
coefficient. In the main component | $9 ), two nucleons 
are combined into spin Snn = and isospin Inn = 1 
so all the two-body subsystems can be in s wave. We 
also allow a mixture of the | $° ) component where both 
spin and isospin are set to be zero, so the orbital angular 
momentum between two nucleons is odd. The J p = 1~ 
state is studied in a similar way as 



J=i ■ 



(aa 1 )- 1 !!^: 



where \$\) denotes S N n = 1 and I N n = 

(Snn — 1 and Inn = 1) component. Note that only the 
main component of | <J)^ -0 ' 1 ) is taken into account in the 
FCA calculation. The wave functions are expanded in 
terms of gaussians in coordinate space, and we minimize 
the total energy of the system with the Hamiltonian given 
below. Detailed explanation of the variational method 
can be found in Ref . [3(| ■ 

We consider the following Hamiltonian in this study: 



H = T + V NN +Re V DN - T c . 



(18) 



where f is the total kinetic energy, Vdn is the DN poten- 
tial term which is the sum of the contributions from two 
nucleons, and T c m . is the energy of the center-of-mass 
motion. For the NN potential Vnn, we use three mod- 
els: HNIR which is constructed from Hasegawa-Nagata 
No.l potential [63], the Minnesota force [Hj], and the 
gaussian-fitted version of the Argonne vl8 potential [65j |. 
The characteristic features of these AW potentials are 
summarized in Appendix [X] For later convenience, we 



where Tn and T c m j jv are nucleonic parts of the kinetic 
term and center-of-mass energy, respectively. 

We take the real part of the DN potential for the en- 
ergy variation, and the imaginary part will be used to 
estimate the mesonic decay width. The energy depen- 
dence of the interaction was treated self-consistently in 
the study of KNN system [30| . While the KN amplitude 
is well calibrated by experimental data such as total cross 
sections and 7rE mass distributions, the DN amplitude 
is only constrained by the mass of the quasi-bound state 
A c (2595) = A*. In addition, the self-consistent treat- 
ment requires some assumption on the energy fraction of 
the DN pair in the three-body system, which cannot be 
determined unambiguously. In this study, therefore, we 
refrain from the self-consistent treatment of the energy of 
the DN subsystem and set the strength of the potential 
at the energy of A* resonance: 



Re v DN (r = 0;W = M A ,) = 



-1336 MeV (I 
-343 MeV (I 



0) 
1)' 
(19) 



with M A . = 2597.1 MeV. In this case, the M A , in / = 
channel is correctly reproduced, while the I — I reso- 
nance disappears, because the strength of the DN po- 
tential (O reduces at the lower energy region as seen in 

Fig.m 

It is useful to introduce one- and two-body densities in 
order to extract the spatial structure of the DNN bound 
state. We first define the one-body densities as 



Pn 



(r)=<*| £ ^(Iri-JlGl-fOm 



p D (r) ={*\8 3 (\r D -RG\-r)\*), 
Pr(r) =PN(r) + Pd(t), 

where Rg is the center-of-mass coordinate of the three- 
body system. The one-body densities represent the prob- 
ability of finding N, D, or any of them at distance r from 
the center of mass of the system. We also define the two- 
body correlation densities as 



Pnn(x) =( * |<5 3 (|ri - r 2 | - x)\ * ), 



i=l,2 



which stand for the probabilities of finding NN or DN 
pair at relative distance x. The root-mean-square radius 



p 3(10(10 
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FIG. 9. Modulus squared of the three-body scattering am- 
plitude for 7 = 1/2 and J = with reduced size of the NN 
radius. 
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FIG. 10. Modulus squared of the three-body scattering am- 
plitude for 1 = 1/2 and J = 1 with reduced size of the NN 
radius 



of particle X, \J (r 2 )x, and relative distance of particles 
X and Y, Rxy, are given as the second moment of the 
one-body and two-body densities, respectively: 

(r 2 ) x = J d 3 r r 2 Px {r), 

r xy = J d 3 x x 2 pxy(x). 

In this setup, since the imaginary part of the DN po- 
tential is not included, the A* appears as a stable bound 
state. Thus, in the variational approach, the DNN 
three-body bound state can be found in the energy region 
below the A*N threshold y/s ~ 3536 MeV. If the three- 
body (quasi-)bound state exists above the A*N thresh- 
old, variational calculation will find the A'*N two-body 
scattering state as the ground state of the three-body 
system. 

A three-body bound state above the ttA c N threshold 
yfs ~ 3363 MeV has a mesonic decay width. The three- 
body decay width can be estimated by the matrix ele- 
ment of the imaginary part of the DN potential as 

where | \& ) is the obtained wave function of the ground 
state. As seen in Fig. [21 the imaginary part of the DN 
potential is much smaller than the real part. This may 
justify the perturbative treatment of the imaginary part, 
which ignores the dispersive effect on the energy of the 
DNN system from the imaginary part. 

V. RESULTS 

A. Quasi-bound states in the FCA approach 

We first study the quasi-bound state found in the FCA 
calculation. In Figs. O and [TU] we show the results for |T| 2 




s" 2 [MeV] 



FIG. 11. Modulus squared of the three-body scattering am- 
plitude for 7=1/2 and J = (with SG) with reduced AW 
radius. 



as functions of the total energy y/s assuming the NN 
system to have reduced size. Both for Inn = 0, Inn — 
1(J = 1, J = 0), we obtain a neat peak. The resonance 
energy for J = is about 3486 MeV and the width is 
extremely small. In the case of J = 1 we have a smaller 
binding and the energy is about 3500 MeV, with a width 
of around 9 MeV. We should note that the binding is 
similar for both the spin channels. The position of the 
peak in this approximation is, in a rough estimate, given 
by the position of the pole of the A c (2595). This gives 
the value of Sj and through Eqs. (|TT|) . (fT!?]) the value of 
s. 

However, one should note the different strength of |T| 2 
in these two cases, but a direct comparison cannot be 
done because the strength of the resonance amplitude at 
the peak is related to the width, which strongly depends 
on the spin. A proper comparison is better done after the 



10 
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FIG. 12. Modulus squared of the three-body scattering am- 
plitude for 7=1/2 and J = 1 (with 8G) with reduced NN 
radius. 



D absorption is included where the widths are similar. 

Next we include the Sg to account for absorption and 
plot \T\ 2 for the DNN system in Figs. QT] and [TJ for 
J = (Inn = 1) and J = 1 (Inn = 0). The difference 
of the peak position by the absorption effect is only a 
few MeV (2-4 MeV) which is certainly within our uncer- 
tainties. The novelty, which is welcome, is that |T| 2 has 
become now wider and acquires a width of about 20-25 
MeV. We are now in a position to compare the strength 
of these two amplitudes and we see that in the case of 
,1 = the strength of \T\ 2 at the peak is about a fac- 
tor 15 larger than that for J = 1. This means that the 
state that we find at J = 1 should be more difficult to 
see, or alternatively we should see the small strength as 
an indication that this state is more uncertain in our ap- 
proximation, as should be the smaller shoulder that one 
can see at higher energies for J = 1 in Fig. [T21 



B. Quasi-bound states in the variational approach 

Now we investigate the same system in the variational 
approach. We first adopt HN1R potential for the nuclear 
force. As a result of the variational calculation, we have 
found that the total spin J = 1 system (Inn = 0) is 
unbound with respect to the A*7V threshold. A bound 
state of spin J = system (Inn = 1) is found at 

B ~ 225 MeV, 

measured from the DNN threshold (~ 3745 MeV). This 
corresponds to the total energy of the three-body system 
as 

M B ~ 3520 MeV. 

We also examine the Minnesota force and Avl8 potential. 
The results are summarized in Table HI together with the 
contributions from the individual terms in Eq. (CG 



TABLE I. Results of the energy compositions in the varia- 
tional calculation for the ground state of the DNN system 
with total isospin 1=1/2 (range parameter a s = 0.4 fm). 
Terms "bound" and "unbound" are defined with respect to 
the A*N threshold. All the numbers are given in MeV. 





HN1R 




Minnesota 


Avl8 




J = 1 


J = 


J = 


J = 




unbound 


bound 


bound 


bound 


JD 


90S 


99^ 


9 r i1 

LU 1 


9DQ 
zuy 


M B 


3537 


3520 


3494 


3536 






26 


38 


22 


-Ekin 


338 


352 


438 


335 


V(NN) 





-2 


19 


-5 


V(DN) 


-546 


-575 


-708 


-540 


T nuc 


113 


126 


162 


117 


Enn 


113 


124 


181 


113 


P(Odd) 


75.0 % 


14.4 % 


7.4 % 


18.9 % 



As seen in the Table U the DNN system in the J = 
channel is bound below the A*N threshold (B ~ 209 
MeV) for all the NN potentials employed^ A large ki- 
netic energy of the deeply bound system is overcome by 
the strong attraction of the DN potential, while the NN 
potential adds a small correction. Comparing the results 
with three different nuclear forces, we find that the bind- 
ing energy is smaller when the NN potential has a harder 
repulsive core (see Appendix [A}. 

In the J = 1 channel, the ground state energy is ob- 
tained slightly above the A*N threshold. The fact that 
the J — 1 channel is unbound is confirmed by changing 
the parameter fi in the trial wave function, which controls 
the size of the total system [30] . By increasing the sys- 
tem size, the total energy gradually approaches the A*N 
threshold. This indicates that the lowest-energy state is 
indeed a two-body scattering state of the A*N channel. 
A large fraction of the odd component in this channel 
(~ 75 %) is realized to enhance the Inn — 1 compo- 
nent which has larger fraction of the Idn = than the 
Inn — component. In fact, pure | (DN)i =n N) state 
can be decomposed into Inn = and Inn = 1 compo- 
nents with the ratio 1:3. Since the Inn — 1 state is the 
odd state in J = 1 (Snn = 1) channel, the 75 % fraction 
of the odd component indicates that the DN pair forms 
the A*. We also examine the J = 1 channel with the 
Minnesota force. Although the repulsive core is soft in 
this case, no bound A*N is found. 

Using the imaginary part of the DN potential, we eval- 
uate the mesonic decay width of the quasi-bound state in 



1 Avl8 case is almost at the A*iV threshold, but we confirm that 
the wave function is localized as we will see in Sec. IV CI 
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TABLE II. Structure of the DNN ground state (range param- 
eter a s — 0.4 fm). \/ (r 2 ) T , y {r 2 ) d and */ (r 2 ) N mean the 
root-mean-square radius of the distribution of total system, 
nucleons and D meson, respectively. -Rjvjv {Rdn) is the mean 
distance between two nucleons (D meson and a nucleon) in 
the DNN. Rdn (I) is the mean distance of a DN component 
with isospin I. All the numbers are given in fm. 







HN1R 

J = 1 


J = 


Minnesota 
J = 


Avl8 
J = 


y/{r 2 )T 




4.81 


0.75 


0.50 


1.26 


^J(r 2 )N 




5.61 


0.88 


0.59 


1.47 


^{r 2 )D 




2.52 


0.41 


0.28 


0.67 


Rnn 




10.04 


1.55 


1.03 


2.62 


Rdn 




7.11 


1.12 


0.76 


1.87 


Rdn(I = 


0) 


4.52 


0.83 


0.62 


1.28 


Rdn(I = 


1) 


10.03 


1.57 


1.03 


2.65 



the J — channel, T^y c n- The results are 20-40 MeV as 
shown in Table fl] This corresponds to the result of FCA 
without the D absorption, where the width is less than 
10 MeV. Note, however, that in the variational approach 
we have evaluated the width perturbatively, while in the 
FCA the evaluation is done nonperturbatively. In this 
sense, T„y c n obtained in the variational approach can 
only be regarded as an estimation of the mesonic decay 
width. 



C. Structure of the DNN quasi-bound state 

To further investigate the structure of the DNN sys- 
tems, we calculate the expectation values of various dis- 
tances of the obtained wave function. The results of 
the root-mean-square radii and the relative distances are 
shown in Table [TTJ Except for the Avl8 case where the 
wave function spreads due to the weaker binding, the 
size of the DNN bound state in the J = channel is 
smaller than the KNN system, in which the NN and 
KN distances are Rnn ~ 2.2 fm and R^ N ~ 1-9 fm. It 
is, on the other hand, acceptable to use the reduced size 
of Eq. dTUJ) for the NN distribution in the FCA calcula- 
tion, given the uncertainty that arises from the choice of 
the NN interaction. The large relative distances in the 
J = 1 channel also reflect the nature of the scattering 
state in this channel. 

In Fig. [T31 we show the one-body densities of the nu- 
cleon and D meson of the quasi-bound state with the 
HN1R potential. It is clear that the D meson distributes 
more compactly than the nucleons. This result indicates 
a schematic picture where the D meson sits at the center 
and nucleons circulates around it. 

It is instructive to look at the DN correlation in more 




r [fm] 

FIG. 13. (Color online) Top: One-body densities pjv(r) and 
Pd(t) in the J = channel with HN1R potential. Bottom: 
the same plot of the densities multiplied by r 2 . 



detail. In Fig. Q31 we show the DN two-body correlation 
density as well as its isospin decomposition. It is seen 
that the 1 = component distributes more compactly 
than the 1=1 component, which reflects the strength of 
the attraction in each channel [see Eq. I|19jl ]. Moreover, 
the 1 = component is similar to the distribution of 
the relative distance of the DN two-body bound state 
PA*(r). This indicates that the structure of the A* is 
maintained even in the three-body system. This feature 
has also been found in the KNN system [3fJ, [66| . 

As in the case of the KNN system, the survival of the 
A* in the three-body system opens the possibility of the 
"A*-hypernuclei", in which the A* is treated as an effec- 
tive degrees of freedom [6^] . In fact, this picture is 
more suitable in the charm sector, since the width of the 
A* is smaller than the A* so the effect of the imaginary 
part in the calculation should be smaller. Note also that 
the binding of the DN system is as large as 200 MeV, 
while the binding of the A* N is much smaller, especially 
for the case of the realistic Avl8 potential. 

We have examined theoretical uncertainties in the con- 
struction of the potential. The range parameter of the 
DN potential a s is introduced in Eq. §5§ and chosen to be 
0.4 fm. When we adopt a s = 0.35 fm, the binding energy 
changes by a few MeV, and the size changes less than 0.1 
fm. The Minnesota potential has a parameter u which 
controls the strength of the NN odd force [64| . The ef- 
fect of the slight inclusion of the odd force (u = 0.95) 
turns out to be very small, less than 1 MeV. We thus 
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r [fm] 

FIG. 14. (Color online) Top: Normalized DN two-body corre- 
lation density pdn{t) with isospin decomposition. The 1 = 
DN bound state (A*) correlation density is also shown for 
comparison. Bottom: the same plot of the densities multi- 
plied by r 2 . 



conclude that these uncertainties are much smaller than 
the dependence on the choice of the NN potential. The 
variation of the values in Tables HI and ITT1 can be regarded 
as the theoretical uncertainties in the present calculation. 



VI. DISCUSSION 
A. Comparison of two approaches 

We have presented the results of two approaches, the 
Faddeev FCA calculation and the variational calculation. 
In the total spin J = channel, both approaches find 
a quasi-bound state around 3500 MeV which is below 
the A*N threshold. The assumed NN distribution in 
the FCA turns out to be similar with that found in the 
variational calculation by minimizing the total energy. It 
is therefore reasonable to conclude that these approaches 
find the same quasi-bound state. 

The spin J = 1 channel, on the other hand, has dif- 
ferences in the two approaches. The lowest-energy state 
obtained in the variational calculation is a A*A^ scat- 
tering state, while a narrow peak is found in the FCA 
amplitude below the A*N threshold, although the sig- 
nal strength is not so significant as the J = case. A 
major reason of this discrepancy may be traced back to 



the DN interaction in the isospin 1=1 channel. In the 
original coupled-channel amplitude, there is an / = 1 
quasi-bound state, which induces the bound state in the 
FCA. As discussed in Sec. Ill B[ however, the energy de- 
pendence of the DN potential in the variational approach 
is fixed at the energy of the A* in the 1 = channel. This 
reduces the strength of the 1 = 1 amplitude, and the two- 
body quasi-bound state is not generated in the effective 
potential. Since the total spin J = 1 channel has larger 
fraction of the 1=1 DN amplitude, this difference is 
enhanced and results in different three-body results. 

In fact, we may artificially adjust the condition l|19p to 
generate a quasi-bound state in the 1 = 1 channel in the 
variational approach. By setting the strength of the DN 
interaction at W ~ 2766 MeV in the 1=1 channel, a 
quasi-bound state is generated in the 1=1 DN channel. 
In this case, the energy dependence of the DN interac- 
tion is fixed at each isospin channel, and the strength of 
the DN attraction is increased in the 1=1 channel. By 
performing the three-body calculation, we find that the 
binding energy in the J = quasi-bound state are in- 
creased by 10-50 MeV, depending on the NN interaction 
employed. This is because of the increase of the attrac- 
tion, and the binding energy appears to be closer to the 
FCA result. In the 3=1 sector, only the Minnesota 
potential supports a bound state with B = 214 MeV, 
while no state is found below the A*N threshold with 
the other two AW interactions. Given the uncertainty in 
the choice of the NN interaction, the present result does 
not strongly support the existence of the quasi-bound 
state in the J = 1 sector. In order to pin down the 
J = 1 quasi-bound state, it is necessary to accumulate 
the experimental information of the DN 1 = 1 scattering 
amplitude, or the information on the negative parity X* 
resonance. 

In addition, we should also remember that the two ap- 
proaches employ different approximations. In the FCA, 
the dynamics of the nucleons is not solved explicitly, 
while the imaginary part of the DN potential is not taken 
into account in the variational approach. In both cases, 
explicit nY c N dynamics is approximated at different lev- 
els (see the discussion in Ref. [36j]), whereas its impor- 
tance has been pointed out in the strangeness sector [3 11 ] . 
These effects can also be responsible for the difference of 
the results in the two approaches. 



B. Comparison with KNN results 

It is instructive to compare the DNN quasi-bound 
state with the corresponding KNN state in Ref. [3(| • In 
both cases, we obtain a quasi-bound state, but the DNN 
system has a larger binding energy and a narrower width. 
This is in parallel with the properties of the DN and KN 
two-body quasi-bound states, and they are closely related 
through the DN and KN interactions. 

As discussed in Sec. Ill Al the D meson can be more 
strongly bound in a nucleus than K meson by two rea- 
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TABLE III. Binding energies of the three-body bound state in 
J — channel measured from the three-body threshold with 
different meson mass and different meson-nucleon potential. 







m = 


m — mo 


V = 




~ 30 MeV 


~ 190 MeV 


V = 


V D 


~ 40 MeV 


~ 230 MeV 



sons. On one hand, the coupling itself is stronger, and on 
the other hand, the heavier mass of the D meson is ad- 
vantageous to increase the binding. So, we can consider 
two hypothetical variants between the DNN system 
(B ~ 230 MeV) and KNN system (B ~ 30_MeVj3; case 
I: kinematics of the DNN system with the KN potential 
(to = ttl d , V = Vjq ), and case II: kinematics of the KNN 
system with the DN potential (m = m^,V — Vr>)- The 
result of the variational calculation shows that B ~ 40 
MeV for case I and B ~ 190 MeV for case II. As summa- 
rized in Table IIII1 the suppression of the kinetic energy 
by the heavy D mass is more important for the strong 
binding of the DNN system. One should note that in 
the present case, the strength of the two-body interac- 
tion is fixed at the energy of the two-body quasi-bound 
state. Since the DN two-body bound state locates 200 
MeV below the DN threshold, the strength of the poten- 
tial is reduced, as seen in Fig. [TJ Thus, in the present 
prescription, the attractive strength of the DN potential 
is not very much different from the KN one, and the re- 
sult of case I does not very much deviate from the KNN 
quasi-bound state. 

The narrow width of the DNN system is a consequence 
of the narrow width of the A* (2595). This is partly be- 
cause of the small transition coupling which is suppressed 
by the exchange of the heavy flavor, but the main reason 
is the suppression of the phase space due to the large 
binding energy. In this sense, the heaviness of the D 
meson is essential to realize the deep and narrow DNN 
quasi-bound state. 

C. Possible experiments to produce the DNN state 

The very narrow width of the DNN system is quali- 
tatively different to the KNN one where the width was 
so large as to make its experimental observation unfea- 
sible. In the present case there is a clear situation and 
there are no problems in principle for the observation of 
the state. In the FCA calculation, we observe that the 
two-nucleon absorption width is larger than the three- 
body decay width. This indicates that the DNN quasi- 
bound state can be more easily seen in the two-baryon 



2 Here we also set the strength of the KN potential at the energy 
of the A* for comparison with the DNN calculation. 



final states such as A c iV. The findings of the present 
work should stimulate efforts to find suitable reactions 
where this state could be found. 

As a suggestion in this direction we can think of the 
p 3 He -> D°D°pn -> D°[DNN] reaction, which could 
be done by FAIR at GSI. With a p beam of 15 GeV/c 
there is plenty of energy available for this reaction and 
the momentum mismatch of the D° with the spectator 
nucleons of the 3 He can be of the order of 550 MeV/c, 
equivalent to an energy of 80 MeV for the D, small com- 
pared with the scale of the binding (> 200 MeV). With 
an estimate of a ~ 10 — 20 nb for pp — > D°D° produc- 
tion [69|, [7(3] one would expect several thousand events 
per day for the background of the proposed reaction [7l[ . 
A narrow peak could be visible on top of this background 
corresponding to the DNN bound state formation. 

Another possibility is the high-energy tt induced re- 
action. An analogous reaction is Tr~d — > D~ D + np — > 
D~[DNN] where the relevant elementary process is 
ir-N -4 D+D-N. Since the DN pair in the DNN 
system is strongly clustering as the A*, the reaction 
ir~d — >• D~A*n — s- D~[DNN] is also another candi- 
date. The elementary reaction ir~p D~A* is relevant 
in this case. Such reactions may be realized in the high- 
momentum beamline project at J-PARC. 

A different strategy is to look for the formation of 
the quasi-bound state in the heavy ion collisions. It 
has been shown that the hadronic molecular states with 
charm quark are abundantly produced at RHIC and 
LHC [zl, |7j| . Although a deeply bound DNN state has 
smaller production yield, it can also be produced via co- 
alescence of the A*N with much smaller binding. A peak 
structure of the DNN state may be seen, for instance, 
in the invariant mass spectrum of the A c 7r _ p or A c p final 
state. 



VII. CONCLUSIONS 

We have studied the DNN system with 1 = 1/2 and 
have found that the system is bound and rather stable, 
with a width of about 20-40 MeV. We obtained a clear 
signal of the quasi-bound state for the total spin J = 
channel around 3500 MeV. 

We have used two methods for the evaluation of the 
quasi-bound state. The first one used the fixed center 
approximation for the Faddeev equations and the second 
one employs the variational approach with hadronic po- 
tentials in coordinate space. The DN interaction was 
constructed in the field theoretical method with channel 
couplings and a unitary approach dynamically generates 
the A* (2595) resonance as a DN quasi-bound state. 

In both cases, we have found a bound state with an 
energy around 3500 MeV in the J = channel. This cor- 
responds to 250 MeV binding from the DNN threshold. 
The J = 1 channel is more subtle, and the precise DN 
amplitude in the 1 = 1 channel is important for a robust 
prediction in this channel. The mesonic decay width of 
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the quasi-bound state turned out to be less than 40 MeV. 
In addition, the D absorption on two nucleons was eval- 
uated in the FCA formalism using a novel method. Al- 
though the absorption process adds several tens of MeV 
to the width, the total width is still much smaller than 
the binding energy. It is found that the DN pair in / = 
channel in the DNN system resembles the wave function 
of the A c (2595) state in vacuum. Thus, the DNN state 
found here can be interpreted as a quasi-bound state of 
A c (2595) and a nucleon. 

The small width of the DNN quasi-bound state is 
advantageous for the experimental identification. The 
search for the DNN quasi-bound state can be done by 
p induced reaction at FAIR, 7r~ induced reaction at J- 
PARC, and relativistic heavy ion collisions at RHIC and 
LHC. 



VIII. ACKNOWLEDGMENTS 

This work is partly supported by projects FIS2006- 
03438 from the Ministerio de Ciencia e Innovacion 
(Spain), FEDER funds and by the Generalitat Valen- 
ciana in the program Prometeo/2009/090. This research 
is part of the European Community- Research Infrastruc- 
ture Integrating Activity "Study of Strongly Interacting 
Matter" (acronym HadronPhysics2, Grant Agreement n. 
227431) under the Seventh Framework Programme of 
EU. T.H. thanks the support from the Global Center 
of Excellence Program by MEXT, Japan, through the 
Nanoscience and Quantum Physics Project of the Tokyo 
Institute of Technology. This work is partly supported 
by the Grant-in- Aid for Scientific Research from MEXT 
and JSPS (Nos. 24105702 and 24740152). 

Appendix A: NN interactions 

Here we summarize the properties of the NN poten- 
tials used in this study with variational calculation. We 
have examined three kinds of NN interactions. Because 
we work in the isospin symmetric limit, the Coulomb in- 
teraction is not included in all cases. 

Hasegawa-Nagata No.l potential [H| has a three-range 
gaussian form with no odd force. The repulsive core is 
as high as 1 GeV. Because the potential was originally 
introduced for the RGM study, when applied to the two- 
nucleon systems, the attraction is too strong to generate 
a bound state in 1 Sq channel and to overestimate the 
deuteron binding energy. In this study, we have reduced 
the strength of the long-range term (middle-range term) 
by factor 0.25 (0.95) and call it HN1R potential. The 
HN1R potential has no bound state in 1 S'o channel and 
reproduces the NN phase shift data, as shown below. 

Minnesota force [64| is expressed by the sum of two 
gaussians, with relatively soft repulsive core. The pa- 
rameters were chosen so as to reproduce the scattering 
lengths and the effective ranges of the JVJV scattering. 
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FIG. 15. (Color online) Coordinate space NN potentials in 
the So channel. 



As a consequence, the deuteron is bound only with the 
s-wave component, so the tensor force is considered to 
be renormalized in the central part. We set the parame- 
ter u = 1 so that there is no odd force, unless otherwise 
stated. 

Argonne vl8 potential [65( is one of the realistic nuclear 
forces with strong repulsive core. As in Ref. [3(|, we used 
the gaussian-fitted version of the potential with central, 
spin-spin, and I? terms. Since the description of the 
deuteron requires the d-wave mixing which is beyond the 
present model wave function, we only consider the S = 
channel with the Avl8 potential. 

In Fig. [T51 we show the spatial form of the potentials in 
the channel. The phase shifts of the JVAf scattering 
in the 1 Sq channel are shown in Fig. \W\ in comparison 
with experimental data. 



Appendix B: Derivation of the three-body 
amplitude in the charge basis 

In this Appendix, we derive Eqs. ([5]) and §§§ in the 
approach of Ref. [361 ] by applying the following strat- 
egy. We evaluate first the D°pp — > D°pp amplitude 
considering charge exchange processes in the rescatter- 
ing of the D meson. The amplitude will contains total 
isospin /tot = 1/2 and 7 t ot = 3/2. Since we only want 
the /tot = 1/2, we evaluate the scattering amplitude for 
^tot = 3/2 in addition, taking the D + pp — > D + pp ampli- 
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FIG. 16. (Color online) AW phase shifts in the So channel 
calculated by the NN potentials. 



tude, and from a linear combination of the two we obtain 
the Jtot = 1/2 amplitude. This strategy was found most 
practical in Ref. [361 ] . 

For the D°pp — > D°pp we define three components of 
the three-body scattering amplitude; 

a) T p , which are called partition functions, which con- 
tains all diagrams that begin with a D° collision with the 
first proton of the pp system and finish with D°pp, 

b) Tex\ which contains all the diagrams that begin 
with a D + p collision on a np system and finish with D a pp, 
and 

c) Tex i which contains all the diagrams that begin 
with a D + n collision on a np system and finish with 
D°pp. 

These amplitudes fulfill a set of coupled equations 
T p = t p + tpGoTp + t ex Gc,T^) 

rp{p) — Mri rp{n) 
-'-ex — L u 0- I ei 

Tj") = t ex + t ex G T p + t ( n) G T$ 



(Bl) 



where the two-body amplitudes are given as t p = 

tD°p,D°pi tex — tjjOpjj+n, t Q P — tp+ p,D+p, and tg ^ = 

t D+n. D+n- The set of equations (|B1[) are diagrammati- 
cally represented in Fig. [T7] 

By taking into account the phase convention = 
— 1 1/2 , — 1/2) in the isospin basis, we can write all the 
former elementary amplitudes in terms of J = 0, 1 
(t(°), for the DN system, and we find 

t~ = i(t<°>-t< 1J ) 
M - +W 

L — I 



4") = A ( i(o) + t (i)). 

Eliminating T$ and T^ in Eq. CEE) we obtain 



(B2) 
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(1 - t p G )(l - t^G t p G ) - tl x t ( g'Gl 
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FIG. 17. Diagrammatic representations of the partition func- 
tions for the D°pp — > D°pp. 



which in isospin basis can be simplified to 

l( f (0) +t (D)_ f (Q) t (i) a G .g 



(i - G t( i ))(i + |(tw - m)G Q - cimtw) 



The total D°pp D°pp amplitude would be 2T p ac- 
counting for the first interaction of the D° with either of 
the protons. 

Now we take into account that in the basis of 

\Itot, l3,tot) 

\D°pp) = -(-^13/2, 1/2) + y||l/2, 1/2)) (B4) 
and thus 

(l/2|T|l/2) = ^({D°pp\T\D°pp) - |<3/2]Tl3/2» (B5) 

The (3/2|T|3/2) amplitude is particularly easy to obtain. 
In this case we take the D + pp — > D + pp transition and 
diagrammatically we have the mechanism of Fig. [TS] for 

(3/2) 

the only partition function Tp . 
Hence 



T^=t^ +t^G TWV 



(B6) 



and the total T^ 3 / 2 ) amplitude will be 2Tp 3 ^ 2 \ accounting 
for the D + interacting first also with the second nucleon. 
We have 



T (s/2) = t = 

l-Go4 p) l-Go*W 



(B7) 
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FIG. 18. Diagrammatic representation of the partition func- 
tion for 1=3/2. 



We can now use Eq. JB5J) and find for the total amplitude 
(including the factor two for first interaction with cither 
proton) 



T (l/2) _ 3T 



1 - G t (1) 



^(0) + IfCD _ I t (D(t(i) - t(°))G - 2t(°hW 2 Gt 
(1 - G i (1) )(l + U t(1) ~ t{0) )°o ~ G§t (0) * (1) ) 



which can be simplified dividing the numerator by (1 
Got^ 1 ') with the final result 



y,(l/2) = 



|t<°> + ±i« + 2G t<°>f< 1 > 
I ~t(°))G -G§i(W 



This corresponds to Eq. (JSJ). 

The case of Sjvat = 1 (Ijvjy — 0); that we also study 
here, can be done in a similar way, but this was done in 
[48| for K~d at rest and in [35| for the K~d interaction 
below threshold. We quote here the formula that was 
obtained in 35] which we use here too. 



T D o 



1_ i( t (i) -t(o))Go-Ggi(W 



This corresponds to Eq. . 
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